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ABSTRACT 

We present a detailed analysis comparing the velocity fields in molecular clouds and the atomic 
gas that surrounds them in order to address the origin of the gradients. To that end, we present 
first-moment intensity-weighted velocity maps of the molecular clouds and surrounding atomic gas. 
The maps are made from high-resolution ^'^CO observations and 21-cm observations from the Lei- 
den/Argentine/Bonn Galactic H I Survey. We find that (i) the atomic gas associated with each molec- 
ular cloud has a substantial velocity gradient — ranging within 0.02 to 0.07 km s^^ pc^^ — whether or 
not the molecular cloud itself has a substantial linear gradient, (ii) If the gradients in the molecular 
and atomic gas were due to rotation, this would imply that the molecular clouds have less specific 
angular momentum than the surrounding H I by a factor of 1 - 6. (iii) Most importantly, the velocity 
gradient position angles in the molecular and atomic gas are generally widely separated — by as much 
as 130° in the case of the Rosette Molecular Cloud. This result argues against the hypothesis that 
molecular clouds formed by simple top-down collapse from atomic gas. 

Subject headings: ISM: clouds — ISM: kinematics and dynamics — ISM: individual (Perseus molec- 
ular cloud, Orion A, NGC 2264, Monoceros R2, Rosette molecular cloud) ~ ISM: 
molecules 



1. INTRODUCTION 

Giant molecular clouds, both Galactic and extragalac- 
tic, are observed to have velocity gradients that many 
authors have interpreted as being caused by rotation 
(e.g., Kutner et al. 1977, Phillips 1999; Rosolowsky et 
al. 2003). If we start with the premise that these clouds 
are rotating because they have inherited the angular mo- 
mentum of the rotating galactic disk out of which they 
formed, conservation of angular momentum should pro- 
vide clues that give us insight to the origin of their for- 
mation. However, simple formation theories that assume 
giant molecular clouds (GMCs) form by condensing out 
of the Galactic disk are at odds with some of the ob- 
servations. For instance, they do not adequately explain 
why the directions of CMC velocity gradients are not 
typically aligned with the direction of Galactic rotation, 
counter to the expectation of conservation of angular mo- 
mentum. Furthermore, simple formation scenarios tend 
to overpredict the observed specific angular momentum o 
f GMCs (e.g.. Blitz 1990; Rosolowsky et al. 2003). This 
is the so-called "angular momentum problem." 

Provided there is no transfer of angular momentum, 
the angular momentum of a CMC should be equal to 
that of the gas out of which it formed. But Blitz (1990), 
working under the assumption that the velocity gradi- 
ents in molecular clouds are due to rotation, showed that 
the angular momentum due to Galactic differential ro- 
tation in the solar neighborhood interstellar medium is 
consistently greater than that contained within molecu- 
lar clouds. Even molecular clouds with the largest ob- 
served velocity gradients — such as the Rosette and Orion 
A molecular clouds — have less specific angular momen- 
tum compared to the ISM from which they presumably 
formed. Blitz (1990) also pointed out that because the 
molecular clouds in his sample are rotating in a sense op- 
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posite to that of Galactic rotation, they could not have 
conserved angular momentum from the initial states he 
calculated, unless the local Galactic rotation curve is 
falling, or unless the clouds always collapsed azimuthall 

y- 

The angular momentum problem also extends to ex- 
tragalactic molecular clouds. Roso-lowsky et al. (2003) 
showed that simple CMC formation theories consistently 
overestimate the magnitude of the observed angular mo- 
mentum of molecular clouds in the galaxy M33. On av- 
erage, they found that simple theory overpredicts the 
observed magnitudes of specific angular momentum by 
more than a factor of 5. Furthermore, they found that 
40% of the GMCs in M33 are counterrotating with re- 
spect to the galactic plane. 

In this paper, we shed light on the angular momen- 
tum problem by doing a detailed analysis of the kinemat- 
ics in GMCs and the surrounding ISM. Whereas previ- 
ous studies estimated the initial angular momentum im- 
parted to GMCs from the Galactic rotation curve (e.g.. 
Blitz 1990; Blitz 1993), we compare the velocity fields of 
GMCs to those of the ISM with which they are associ- 
ated directly from observation. In light of the observa- 
tion that GMCs have a spatial and kinematic correlation 
with high-surface density atomic gas (§3), we pose the 
question: Does the rotation of the large-scale H I associ- 
ated with GMCs mirror that of the GMCs themselves? 
Our primary goal is to determine whether or not rotation 
is the cause of the velocity gradients in GMCs. To that 
end, we create first-moment maps of the molecular clouds 
in our sample and of the atomic gas surrounding them 
are for comparison. In the following section, we describe 
the ^'^CO and 21-cm data used to conduct this study. In 
§3, we create intensity-weighted first moment ^'^CO and 
H I maps of five Galactic clouds: Perseus, Orion A, NGC 
2264, Monoceros R2 (MonR2), and the Rosette. The re- 
sults from these measurements are given in §4, and a 
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summary is provided in §5. 

2. DATA 

To measure velocity gradients and other properties 
across the molecular clouds, we use high-resolution, 
high-sensitivity published ^^CO observations. Because 
'^'^CO emission is nearly always optically thin in Galac- 
tic GMCs, we have the advantage of getting a detailed 
view of the kinematic structure of the molecular clouds in 
our sample. And because ^^CO has narrower line widths 
than the optically thick ^^CO, the former permits finer 
separation of velocity components than the latter. Nev- 
ertheless, the large scale velocity gradient of a GMC mea- 
sured using ^^CO is generally consistent with that mea- 
sured using ^^CO, since the gradient is being measured 
across several parsecs and small scale variations in the 
velocity field tend to get averaged out. 

The data for Perseus, Orion A, NGC 2264, and the 
MonR2 molecular clouds were generously provided by 
J. Bally (see Bally et al. 1987). Observations of these 
clouds were taken at the AT&T Bell Laboratories 7 m 
telescope and have a beam size of 100". Perseus, Orion 
A, and NGC 2264 data were resampled onto 60" grids. 
The MonR2 data were resampled onto a 30" grid. The 
spectral resolution of 128 channels at 100 kHz corre- 
sponds to a velocity resolution of 0.27 km s~^. On the T4 
scale, the cubes have rms noise levels of 0.17 K (Perseus), 
0.32 K (Orion A), 0.42 (NGC 2264) and 0.29 K (MonR2). 

J. Williams and M. Heyer graciously provided the 
FCRAO data of the Rosette Molecular Cloud (see Heyer, 
Williams, & Brunt 2006). The beam size is 47" and the 
data were interpolated onto a 20" grid. The spectral res- 
olution is 59 kHz per channel, and the velocity resolution 
is 0.133 km s~^. The data have an rms noise of 0.21 K 
in T^, similar to that of the Bell Labs data. 

The H I data are obtained from the Lei- 
den/ Argentine/Bonn (LAB) Galactic H I Survey 
(Kalbcrla ct al. 2005), which spans velocities from 
-400 km s-i to -)-400 km s'^ (The LSR velocities of 
the molecular clouds in our sample range from 6 to 14 
km with velocity dispersions of a few km s"-'^.) The 
survey has a half power beam width of 0.6°, velocity 
resolution of 1.3 km s^^, and an rms noise level of 0.07 
K. The high sensitivity and resolution of the LAB data 
set enables a detailed study of the atomic gas from 
which the GMCs formed. 

3. ANALYSIS 

Observations in the Milky Way indicate the molecular 
clouds are typically associated with high-density atomic 
gas with column densities around iV(H I) ~ 2 x lO^^cm"^ 
(e.g., McKee & Ostriker 2007). The Rosette Molecular 
Cloud is a prototypical example: Williams ct al. (1995) 
measured the mean column of the H I associated with 
the GMC to be 1.3 x lO^^cm^^. Furthermore, in ex- 
ternal galaxies, GMCs are often observed to be located 
on or near bright H I peaks. In the Large Magellanic 
Cloud, Mizuno et al. (2001) observed that most GMCs 
are associated with H I having column densities greater 
than 10^1 cm"2 Rosolowsky et al. (2003) found that 
every GMC they identified in M33 lies on an over-dense 
H I filament, though every over-dense H I region does 
not contain a GMC. This seems to imply that high col- 
umn density atomic gas is necessary but not sufficient for 



GMC formation. 

It is our goal to do a detailed comparison of the velocity 
fields in the GMCs and local atomic gas. To that end, we 
first describe how the physical properties — including col- 
umn densities, masses, and velocity gradients — in both 
the atomic and molecular gas are determined. We then 
provide our criteria for choosing spatial and kinematic re- 
gions of atomic gas associated with the molecular clouds. 
Lastly, in this section, we discuss how we estimate the 
specific angular momentum. These parameters are sum- 
marized in Tables 1 and 2. 

3.1. Cloud Properties 

The H I column density, A''(H I), is calculated along 

each line of sight by integrating the atomic hydrogen 
emission above a certain background value and over the 
selected velocity range, (see §3.1 and Table 1), 



iV(H I) = 1.82 X 10^^ 



b.HI 



K km ^ 



— dv cm ^, (1) 



where Tb.Hi is the brightness temperature of the H I ob- 
servations and dv is the channel velocity width. Assum- 
ing that the H I is optically thin. Equation 1 provides a 
lower limit to measured column density, which we con- 
vert into units of surface density in Figures 1 5. 

The total H I mass is then determined by summing 
over all pixels in the map where emission is detected, 
that is, where the emission is at least three times the 
root-mean-square (rms) noise level: 

Mhi = ^ M • mn • iV(H I)pixei • {(fAaA6), (2) 

pixels 

where /z = 1.36 is the correction for helium, d is the 
distance to a given molecular cloud, and d'^AaAS is the 
area of one pixel, which corresponds to one resolution 
element. 

The i^CO column density, iV("CO), is derived assum- 
ing that the ^^CO emission is optically thin and in local 
thermodynamic equilibrium. Following Frerking et al. 
(1982), 



Ar(i3C0) = 2.13 X 10" [1 - e-5-287/Tex]-i 



Th,co dv 
K km s-i 



cm 



(3) 



where / Tb.co dw is the integrated ""^^CO intensity and 
Tex is the excitation temperature. Normally, Tex is de- 
termined by measuring the ^^CO radiation temperature 
toward ^^CO peaks. Since we lack ^^CO observations at 
the same resolution as the ^^CO data, we use a constant 
excitation temperature in the calculation of (-'^^CO) for 
each of the five GMCs. Based on the following argu- 
ments, we adopt a value of 20 K for each of the GMCs. 
If the actual excitation temperature in a given region is 
between 10 and 30 K, the derived ^^CO column density 
will be in error by less than a factor of 2. 

Castets et al. (1990) showed that the ^^CO emission in 
Orion A mainly arises from regions where Tex ~ 20 — 25 
K, and dense cloud cores have temperatures of Tox ~ 
15 — 20 K. Nagahama et al. (1998) showed that, with 
the exception of two peaks at I ~ 209° and I ~ 212.5° 
associated with embedded young stellar groups, Tex rises 
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slowly and monotonically in Galactic longitude, with an 
average ranging from 13 — 20 K. In his study of molcciilar 
clouds, including Perseus, Orion A, and MonR2, Carpen- 
ter (2000) adopts a constant value of Tex = 10 K, though 
the coefficient in his formula for A^(^''CO) yields slightly 
higher values than ours in Equation 3. Thus, for Orion A, 
Perseus, NGC 2264, and Monoceros R2, we use Tgx = 20 
K in our calculation of iV(^^CO). 

Williams et al. (1995) showed that Tex in the 
Rosette decreases slowly with increasing distance from 
the Rosette Nebula (centered at I = 206.2°, b = -2.1°) 
from ~ 20 to 5 K. In our analysis, we adopt a uniform 
value of Tox = 20 K for the Rosette. Our estimate of the 
cloud's mass using this value is slightly lower than that 
estimated by Williams et al. (1995), who measured the 
mass over a larger surface area (see below). 

Next, the H2 column density is evahiatcd assuming 
a ratio of N{}i2)/N{^^C0) = 7 x 10'^ (Frerking et al. 
1982). Pixels having values at least three times the rms 
noise level are counted as detected. Finally, the molecu- 
lar mass M13 is calculated over the areas where emission 
is detected (that is, higher than the 3-crrms level) using 



M 



13 = J2 1^' "^H2 ■ A^(H2)pixei ■ (d'^AaAS), (4) 

pixels 

where is the mass of an H2 molecule. We note that 
the areas over which emission is detected and, subse- 
quently, the masses we calculate will be smaller than 
cited in previous studies in which these quantities were 
measured using '^^CO emission. This is because the 
stronger -'^^CO line is observed over larger areas in GMCs 
than the ^^CO line. In the Rosette, for instance, we mea- 
sure a projected area of 1500 pc^ and a mass of 6.0 x 10^ 
Mq, while WilUams et al. (1995) measure 2200 pc^ and 
7.7 X 10" Mq. In Table 1, both the ^^CO masses and pre- 
viously measured ^^CO masses are listed for the GMCs. 

3.2. Velocity Gradients 

Velocity gradients are measured from first moment 
maps of the atomic and molecular gas. First, the 
intensity-weighted average velocity along each line of 
sight is determined using 



(5) 



where Vi and Ti arc the velocity and temperature at loca- 
tion i. Following Goodman et al. (1993), the uncertainty 
of a given measurement is 
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1.2 



' peak 



{dv Aiifwhm) 



1/2 



(6) 



where Trms is the spectrum noise, Tpeak is the maximum 
temperature along the line of sight, and Awfwhm is the 
FWHM linewidth of the spectrum along the line of sight. 

A plane is then fitted to the first moment map of ve- 
locity centroids, as in Goodman et al. (1993), assuming 
a linear velocity gradient: 



Vis 



Vo + a{x - xo) -)- h{y - yo), 



(7) 



where is the mean cloud velocity, and (xq, J/o) is an ar- 
bitrary reference position, which we take to be the center 



of our maps, and the coefficients are 

dv ^ dv 
dx' dy' 

The gradient magnitude and direction, fl and 0, are de- 
rived from the coefficients to the fit. 
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where 6, measured in degrees East from North, points in 
the direction of increasing velocity. Note that since we 
have no information regarding the inclination of a given 
cloud, i, to our line of sight, gradient measurements are 
underestimates of the actual values, i^true = sini. 

The imcertainties in these values are calculated by 
propagating the errors in the coefficients. To check 
whether planes are good fits to the velocity centroid 
maps, we make plots of the central velocity at a given 
location in the cloud versus the perpendicular offset from 
the cloud's rotation axis. This is done by taking the av- 
erage velocity along lines parallel to the rotation axis at 
various distances. In most cases, as will be discussed 
in §4, these plots show that planes are good fits to the 
velocity centroid maps. 

3.3. Selecting H I Regions 

We select H I regions in the position- velocity LAB data 
cube that are centered, spatially and kinematically, on 
the five GMCs in our sample. We do not know a pri- 
ori the kinematic or spatial extent of the H I regions 
that are associated with each GMC nor the extent of 
the H I velocity gradients. Molecular clouds have well- 
defined boundaries at which the molecules are dissociated 
by UV radiation and where there is a distinct transition 
from primarily molecular to primarily atomic gas (e.g. 
Savage et al. 1977; Blitz & Thaddeus 1980). The atomic 
gas associated with GMCs does not have such distinct 
boundaries, however, making it difficult to distinguish 
H I that may be related to GMCs from background emis- 
sion. Thus, we begin by examining regions in position- 
velocity space that are far from the center of the GMC. 
While we want to capture the full extent of any linear 
velocity gradient we may measure in a given region of 
atomic gas, we do not want to make our aperture so large 
that we end up including in our measurements too much 
atomic gas that is unrelated to the molecular clouds. 

We start by varying the size of the region (that is, the 
subcube extracted from the LAB data) and examine how 
the velocity gradient magnitude and direction change. 
We fix the velocity range (see below) and vary the spatial 
size of the region centered on the GMC from about 10 
to several tens of parsecs, in increments of 10 pc. We 
find that the gradient direction remains roughly constant 
until the radius of the region over which it is measured 
reaches 40 ± 10 pc, independent of the size of the GMC. 
Beyond this, the measurements start to fluctuate, as the 
gradient in the vicinity of the molecular cloud becomes 
washed out by unrelated features. 

Thus for each cloud, we end by selecting a spatial 
boundary of atomic gas which extends roughly 30 - 50 pc 



4 



from the center of the CO emission. Because the peaks 
of the H I regions are included in the maps we generate, 
varying the size of the region within 40 ± 10 pc does not 
change the gradient direction by more than a few de- 
grees. The 1-(T uncertainty level of 9 ranges from 3 to 
6 degrees for the clouds in our sample. Our criterion is 
supported by previous studies such as that of Anders- 
son et al. (1991) who found that the spatial extent of 
high-intensity H I halos, measured from the edge of the 
molecular cloud, ranges 5 - 10 pc (see also §5). Figures 
1 through 5 show the H I surface density maps derived 
from the zeroth-moment intensity maps, with the out- 
line of the molecular clouds overlaid. Also overplotted 
in each figure is a dashed circle indicating the H I region 
selected for the analysis. Keep in mind, we are showing 
the ^^CO emission of the molecular clouds, and so the 
maps in Figures 1 - 5 do not show the full extent of the 
molecular emission in the GMCs. 

To choose relevant velocities of the atomic gas asso- 
ciated each GMC, we begin by examining the H I emis- 
sion in the velocity range ±20 km s~^ centered about the 
mean LSR velocity of the ^'^CO emission. Again, this is 
because we want to be sure that our measurements in- 
clude as much as possible of the associated atomic gas. In 
the direction of the Rosette, for instance, the H I has been 
observed to extend several km s^^ beyond the CO emis- 
sion (Williams et al. 1995). Studies of both Milky Way 
molecular clouds (e.g., Wannier et al. 1983, Williams et 
al. 1995) and extragalactic clouds (e.g. Engargiola et 
al. 2003) have shown that the H I emission line profile 
tends to peak in the direction of GMCs. In effect, we are 
using the H I velocity as a proxy for distance in order to 
associate the atomic gas with the GMCs. As Figures 1 
- 5 show, although the H I emission line is broader than 
the CO line, the velocity difference between the peaks in 
the respective lines never exceeds cthI) where cthi is the 
velocity dispersion of the H I profile. 

The bottom panel in Figures 1 through 5 show the 
average ^^CO spectrum through each GMC with the 
H I spectrum in the same direction (within the dashed 
circle) overplotted. Because the '^'^CO-H I peaks arc 
nearly coincident in each case, this indicates that most 
of the H I in the direction of a given cloud is associated 
with that cloud within the selected velocity range. In the 
cases of Perseus, Orion A, and MonR2, the H I emission 
drops abruptly beyond ±15 km s~^ of the ^^CO emis- 
sion. NGC 2264 and the Rosette have more complicated 
H I spectra, each showing double peaks that may be in- 
dicative expansion or of a blended, possibly unrelated 
component. The latter explanation would not be sur- 
prising since, of the five clouds in the sample, NGC 2246 
and the Rosette are located closest to the Galactic plane 
where line-of-sight blending is more of a problem. 

Based on Figures 1 - 5, we determine the boundaries 
of the H I emission we will use for the subsequent analy- 
sis. For Perseus, Orion A, and MonR2, the H I line pro- 
files are approximated as Gaussians and we assume that 
H I emission having velocities within ±2ithi is associated 
with a given molecular cloud. Since both NGC 2264 and 
the Rosette each have a second peak in their H I temper- 
ature profiles (Tb,Hi) at higher velocities and because we 
want to be careful to exclude as much unrelated emis- 
sion as possible, we set a slightly more stringent crite- 



rion on the velocities we select. For each cloud, a maxi- 
imim velocity is identified at the local minimiim in Tb^i 
where the Gaussians overlap. In the NGC 2264 spec- 
trum, for instance, Tb,Hi drops to 29 K at 15 km s~^ and 
then peaks again at around 20 km s""'^ (see Figure 3). 
Thus, we eliminate all emission having velocities above 
15 km s^^, the local ion of the local minimum. All of 
these selections are listed in the third to last column of 
Table 1. 

3.4. Specific Angular Momentum 

Once the magnitudes of the velocity gradients in the 
molecular clouds and the surrounding H I are measured, 
we may calculate and compare their specific angular mo- 
menta, under the assumption that the linear gradients 
are due to solid body rotation. The specific angular mo- 
mentum, j, is simply the total angular momentum of a 
body divided by its mass, 

j = /3fii?2, (11) 

where R is the radius of the region, and the constant 
13 takes into account the moment of inertia of a rotat- 
ing body. For roughly spherical GMCs having constant 
surface mass density distributions, /? = 2/5. Unless oth- 
erwise stated, we take the size of a given molecular cloud 
to be its effective radius, as defined by its projected area: 

-Reff = ^/A/n. 

We would like to estimate the expected specific angu- 
lar momentum initially imparted to a GMC by the ISM 
from which it forms. This depends on the process of 
GMC formation, and we assume here that GMCs form 
via a "top-down" formation mechanism. For instance, 
it has been suggested that molecular cloud formation oc- 
curs when an instability triggers collapse or condensation 
from the Galactic disk (e.g., Mouschovias et al. 1974; 
Blitz & Shu 1980; Elmegreen 1982; Kim et al. 1998). 
Bhtz & Shu (1980) show that "bottom-up" formation 
of GMCs via the random agglomeration of pre-existing 
low-mass clouds is unlikely because of the long timescales 
for this process. Gravitational instabilities and magneto- 
gravitational instabilities, however, tend to proceed more 
quickly. 

If the H I surrounding the GMCs is reflective of the 
ISM out of which the GMCs initially formed, the quan- 
tity (l/2)fiii^ provides an estimate of the initial specific 
angular momentum of the GMCs, assuming they are ro- 
tating. Henceforth, we will often refer to this as the 
expected specific angular momentum. Note that this dis- 
regards the possible effects of magnetic fields; that is, 
we are assuming that the magnetic field strength of the 
forming cloud is 0. The initial angular velocity imparted 
to a given GMC is the local value of $1, which we calculate 
from the first-moment maps of the H I using Eqiiation 9. 
In principle, the size of the region from which a forming 
cloud gathers material, the "accumulation radius" Ra, 
could have a range of values because it depends on details 
of the formation mechanism, as well as on assumptions 
regarding the initial surface density and geometry of the 
gas from which a GMC formed. Nevertheless, we can set 
an effective lower limit on Ra by calculating the size of 
the region from which a GMC must contract to obtain 
its present mass, Mgmc- 

Following Blitz (1993), let us assume that the initial 
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geometry of the collapsing region is a cylinder with a 
diameter equal to its height. The size of the cylinder is 
determined by requiring that the mass contained within 
it is equal to the present mass of the GMC: 

27rpHii?i = Mgmc, (12) 

where pni is the mean volume density of the atomic gas 
from which the GMC formed. Blitz (1993) estimated Ra 
using the mean value of Shi in the Galactic plane near 
the Sun, 5 Mq pc"^ (Henderson et al. 1982). This cor- 
responds to a mass density of pm = 0.0125 Mq pc~^ 
(or a number density of 0.5 cm^'^), using an effective 
scale height of atomic gas of 200 pc in the solar vicin- 
ity (Falgarone & Lequeux 1973). However, using larger 
values of Ehi (leading to lower estimates of Ra) might 
be more appropriate, given the observation that GMCs 
tend to form in regions of H I with densities higher than 
global galactic values (e.g., Engargiola et al. 2003; Imara 
et al. 2010). We estimate the accumulation radii using 
Shi = 10 Mq pc~^, a value more in keeping observations 
of the atomic gas associated with GMCs. For instance, 
around the Rosette, Williams et al. (1995) measure an 
H 1 column density of 1.3 x 10^^ cm^^, which corresponds 
to a Shi ~ 10 Mq pc^^, twice the mean Galactic value. 
From the Sancisi (1974) study of atomic gas near Perseus, 
the inferred surface density is 11 Mq pc~^. Our esti- 
mates of Ra are listed in Table 2, as well as the effective 
radii of the GMCs. 

Finally, we note that the predicted specific angular mo- 
mentum in the atomic gas has a marked dependence on 
the GMC mass. Substituting Ra ~ (Mqmc/phi)^/^ 
from Equation 2.12 for the radius in Equation 2.11 

2/3 

yields jm ~ M^y^^. The masses we calculate using 
the ^^CO observations underestimate the total molecu- 
lar mass of the GMCs. For this reason, we use the larger 
GMC masses, as previously measured using ^^CO obser- 
vations, to predict the initial specific angular momen- 
tum imparted to forming molecular clouds. Typically, 
Ra ^ 3iiGMC- 

4. RESULTS 

Figures 6 through 11 show the H I velocity maps 
(grayscale) overlaid with velocity maps of ^^CO (color) 
for each GMC. Overplotted on these maps are axes of 
rotation: the lines perpendicular to the gradient direc- 
tions, 6qmc and 9hi, as calculated from Equation 10. 
These lines are the position angles of the rotation axes, 
ip, (where t/j = 9 + 90°), of the molecular and atomic gas 
if the gradients are in fact due to rotation. 

In order to check whether planes are good fits to the 
first-moment velocity maps, we make position- velocity 
cuts parallel to the maximum gradient directions and 
plot the results (bottom panels of Figures 6 - 11). We 
are essentially plotting the central velocity at a given 
location in the cloud versus displacement along the gra- 
dient on a pixel by pixel basis. Figures 6, 10, and 11 
show that planes are good fits to the first-moment maps 
of H I surrounding Orion A, MonR2, and the Rosette be- 
cause there is a clear linear trend in the gradient in the 
atomic hydrogen. The velocity fields of the H I associ- 
ated with Perseus and NGC 2264 appear to have more 
complex structure. In the case of the NGC 2264 and 
MonR2 molecular clouds (Figures 9 and 10), however, 



the position-velocity plots do not have monotonically 
increasing or decreasing slopes, indicating a more com- 
plex velocity structure in the molecular gas (red points). 
In these two cases, we nevertheless overplot the lines 
perpendicular to the gradient direction calculated from 
Equation 10. 

The specific angular momenta in the atomic and molec- 
ular gas, listed in Table 3, are compared in Figure 12. 
There appears to be a reasonable correlation — jhi and 
icMC increasing together as jhi oc Jgmc — although 
small number statistics prevent us from making a firm 
conclusion. In each case, the initially expected specific 
angular momentum, jhi, is always greater than jcMC- 
These measurements alone are consistent with a pic- 
ture whereby GMCs form via some top-down mechanism, 
such as a gravitational instability, and somehow shed an- 
gular momentum in the process. But since angular mo- 
mentum is a vector quantity, this scenario is difficult to 
reconcile with the observation that the gradient position 
angles in the molecular and atomic gas differ and appear 
to be uncorrelated (Figure 13). 

The key finding of our analysis is that the regions of 
atomic gas associated with molecular clouds have linear 
velocity gradients, yet the directions of these gradients 
are with one exception — unaligned with the direction 
of the gradients in the associated GMCs. Under the hy- 
pothesis that the gradients are caused by solid body ro- 
tation, this would imply that GMCs are not corotating 
with the surrounding ISM. The second key result is that 
the magnitudes of the velocity gradients in the GMCs are 
larger than the gradient magnitudes in the atomic gas. 
Thirdly, if the gradients in the molecular and atomic gas 
are from rotation, the angular momenta in the molecu- 
lar clouds is less than predicted from calculations of the 
angular momenta in the associated atomic gas. 

Below, we describe the results in detail for each GMC. 
The gradient directions, magnitudes, and specific angular 
momenta of both the molecular and atomic gas are given 
in Table 3. 

Orion A 

Figure 6 shows that the Orion A molecular cloud has a 
large-scale gradient whose velocity decreases from about 
12 km s^^ to 3 km with increasing Galactic longi- 
tude. Of the GMCs in this study, observations of the 
kinematics and morphology of the Orion A molecular 
cloud seem to make the best case for a top-down pic- 
ture of GMC formation. The direction of the gradient 
and the long axis of Orion A are parallel to each other 
and to the Galactic plane. If the gradient were due to 
rotation, this would imply that Orion A is rotating in a 
sense nearly opposite to that of the Galactic disk. This 
result is in agreement with previous studies (e.g., Kut- 
ner et al. 1977, Blitz 1993). The velocity gradient of the 
H I in the immediate surroundings of Orion A points in 
nearly the same direction as that in the molecular cloud. 
Figure 6 shows that the gradient in the atomic hydrogen, 
integrated over the velocity range from about —7 to 22 
km s^^, differs from th e gradient in the molecular cloud 
by only 9°. 

The origin of the velocity gradient in Orion A has been 
debated; it has previously been explained by cloud rota- 
tion (Kutner et al. 1977), expansion driven by the Orion 
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OB association (Bally ct al. 1987), and expansion driven 
by stellar winds from newborn stars (Heyer et al. 1992). 
If due to rotation, the 0.22 km pc~^ gradient in the 
molecular cloud implies a specific angular momentum of 
42.2 pc km s^^. However, considering that Orion A is 
filamentary and much more closely resembles a cylinder 
rotating about its minor axis than a sphere, this esti- 
mate of j is likely to be low limit, since the configura- 
tion of the former has a higher moment of inertia than 
the latter. For a cylinder rotating about its minor axis, 
.7cyi = 1/I2r2(3i?^ + Z/^), where R is the radius of the 
cylinder and L is the length. The cylindrical radius and 
length of Orion A are approximately 1.5° and 9°, corre- 
sponding to 11 pc and 65 pc. Consequently, jcyi = 84 
pc km s~^, much c loser to, but still a bit less than the 
estimated angular momentum of the associated atomic 
gas (107 pc km s~^). 

These two results — that the gradients of the molecular 
and local atomic gas are in near alignment, and the mag- 
nitudes of the specific angular momenta are within range 
of each other — suggest that Orion A may be a case in 
which the molecular cloud and surrounding H I are coro- 
tating and in which conservation of angular momentum 
is demonstrated. As we will see below, it appears to be 
an exceptional case. 

Perseus 

Figure 7 shows that a strong linear gradient exists in 
the Perseus molecular cloud. The velocities in the GMC 

range from about 1 to 11 km s^^ with increasing Galac- 
tic latitude. The gradient direction is tilted about 20° to 
the North- West. Of the molecular clouds in the sam- 
ple, Perseus has the largest velocity gradient at 0.23 
km s~^ pc~^. If the gradient is caused by rotation, from 
Equation 11, the specific angular momentum in Perseus 
is 24.2 pc km s~^, as indicated in Table 3. 

Examination of the black points in the position 
velocity-plot shown in Figure 7 indicates that the veloc- 
ity gradient of the H I region centered on Perseus is not 
linear. As with each cloud in this study, we take the cen- 
ter of the H I field to be located at the position centroid 
of the GMC in Galactic coordinates, (/q, bo). This affects 
the appearance of the position- velocity plot, which shows 
the velocity of points in the field as a function of distance 
from the gradient judged from /qj ^o- The H I surface den- 
sity map in Figure 2 shows a high-surface density (~ 16 
Mq) H I filament extending vertically near a Galactic 
longitude of 157.5°. Whereas the other GMCs in this 
sample have some portion of their molecular material 
laying directly on top of an H I peak, this is not strictly 
the case with Perseus, which is tilted toward the East 
from the bright H I filament. Thus, we generate another 
first-moment map of the H I ce ntered on the filament at 
lo,b() — 157.5,-18, as well as another position-velocity 
plot centered here. Figure 8 shows that a linear fit is a 
much more suitable fit to the gradient in the H I when 
we shift the reference point. When measured this way, 
the magnitude of the gradient measured using Equation 
9 shifts from 0.038 to 0.067 km s^^ pc^^, the value we 
cite in Table 3. The gradient direction is nearly 100°, 
very close to the sense of Galactic rotation, in which ve- 
locities increase from West to East. But the directions 
of the gradients in the molecular cloud and the atomic 



gas surrounding the GMC differ by nearly 120°. 

Of the clouds in this analysis, Perseus has the largest 
velocity gradients in the molecular gas. The magnitude 
of the gradient in the molecular gas (0.23 km s pc~^) 
is 3.4 times that in the H I (0.067 km s^^ pc^^). The es- 
timated specific angular momentum in the atomic gas is 
JHI = 60 pc km s^^, 2.5 times j measured in the molec- 
ular cloud. Although this is consistent with a scenario 
in which the Perseus molecular cloud formed by collaps- 
ing out of the atomic gas, its angular momentum being 
redistributed somehow in the process, based on the ob- 
servation that there is such a large difference between the 
gradient directions in the molecular and atomic gas, it is 
difficult to see how the GMC could have formed in such 
a simple way. 

NGC 2264 

The first-moment maps and position- velocity plot in Fig- 
ure 9 show that the NGC 2264 molecular cloud and 
surrounding H I both have complex kinematic features. 
Though there is no significant linear gradient across the 
entire GMC, there is a weak gradient in the atomic gas. 
The gradient in the atomic gas is stronger to the South 
of the rotation axis, where the velocity decreases from 
about 8 to 5 km s~^ as the Galactic latitude increases 

from 1.8° to +1.5°. As seen in Figure 3, the NGC 

2264 molecular cloud is composed of two main structures. 
The structure in the North has a larger range of veloci- 
ties, with V = 2 - 10 km s~^, than does the structure in 
the South, which averages around 6 km s~^, close to the 
mean velocity of the molecular cloud as a whole. This is 
also the average velocity of the H I peak located at near 
I = 201.5°, b = 0.5° (Figure 3), which is associated with 
the southern segment of t he GMC. 

The magnitude of the gradient in the atomic gas near 
NGC 2264 is 0.019 km s^^ pc^^ which is less than 
fl{= 0.025 km s^^ pc^^) in the solar vicinity assuming 
a fiat rotation curve. The magnitude of the gradient fit- 
ted across the entire GMC, using Equation 9, is 0.046 
km s~^ pc'^. But the northern segment of the cloud in 
isolation has a gradient closer to 0.08 km s~^ pc~^. 

It seems unlikely that rotation is the origin of the ve- 
locity field in the NGC 2264 molecular cloud. Table 3 
records the velocity magnitude and direction calculated 
by fitting a plane to the field, even though the gradient 
is not linear. Using these values, the specific angular 
momentum of the GMC, 4 pc km s"-'^, is smaller than 
the expected value by at least a factor of 6. The gradi- 
ent direction in the GMC, pointing nearly 50° East from 
North, differs from 9m by 140°. An alternative explana- 
tion of the velocity field in NGC 2264 may be, at least in 
part, the internal stellar activity. It is noteworthy that 
the H II region associated with NGC 2264 (the Cone Neb- 
ula) is located in the north, near the part of the GMC 
that has a larger velocity dispersion than the southern 
part. The nebula, at I = 202.95, b = 2.20 (Kharchenko 
et al. 2005) is located between two high velocity regions 
of the GMC that are receding at speeds near 10 km s~^. 
This morpholog y is suggestive of an expanding ring seen 
from an edge-on perspective. Taken all together, it ap- 
pears that these features could be causally connected: an 
H II region causes the high-speed expansion of the sur- 
rounding molecular gas, which subsequently sweeps up 
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atomic gas into a high-density ridge. It is perhaps sig- 
nificant that we see a similar pattern in the Rosette. 

Yet another explanation of the kinematics of NGC 2264 
is put forth by Furesz et al. (2006), based on their find- 
ing that the stars and ^''CO emission in NGC 2264 are 
well-correlated in position-velocity space. They suggest 
that the velocity field of the GMC is explained by the 
models of Burkert & Hartmann (2004), in which molec- 
ular clouds form from supersonic collisions of gas. This 
argument seems to be corroborated by the overall pat- 
tern of the cloud's velocity structure: higher velocities at 
the outskirts of the cloud and lower velocities proceeding 
toward the center of the cloud. This pattern could be 
explained by gravitationally driven infall at the interface 
of the colliding flows. In this scenario, star formation oc- 
curs preferentially in the condensations that develop dur- 
ing the collapse of the molecular cloud, thus explaining 
the position-velocity correlation observed between stars 
and high-density molecular gas. 

The colliding flows hypothesis may also explain the 
velocity structure of the atomic gas. Figure 9 shows 
that the H I has slightly higher velocities at positions 
in the field far from the gradient axis and lower veloci- 
ties close to the axis — again, suggestive of infall acceler- 
ated by gravity. If other molecular clouds form by this 
mechanism, however, it is unclear why we do not see this 
pattern in the atomic gas surrounding the rest of the 
GMCs in our sample. Possibly, the NGC 2264 system 
still retains signs of its early formation history due to 
the relatively young age of the GMC (^1 3 Myr; Flac- 
comio et al.' 2000; Ramirez et al. 2004). It is' difficuh 
to make a definitive conclusion since we are observing 
H I so close to the Galactic plane where the problem of 
line-of-sight blending is exacerbated. 

MonR2 

Figures 4 and 10 show that the MonR2 molecu- 
lar cloud sits right on top of a high-density cloud of 
H I that has a strong linear velocity gradient. The gra- 
dient, which extends for over 80 pc. has a magnitude of 
0.035 km s^^ and increases from about 7 km s~^ to 12 
km s^^ from South to North- West. If the gradient is due 
to rotation, j = 47 - 180 pc km s"^, depending on the 
size of the region considered. From conservation of angu- 
lar momentum, one might expect MonR2 to have a gra- 
dient somewhere between SI = 0.38 - 0.84 km s~^ pc~^, 
pointed about 11° West of North like the gradient in the 
atomic gas. This is not what we see in Figure 10, how- 
ever. 

The velocity field of the MonR2 molecular cloud has a 
much more complex structure than does the atomic gas. 
The outer northern, western, and southern edges of the 
cloud reach velocities up to 12 13 km s~^. The inner 
portion of the cloud is moving at slower speeds around 
9-11 km s^^, which is also the velocity range at which 
the H I peak located at I = 213.5°, b '= -12.5° domi- 
nates. Considering the apparently random nature of the 
velocity field, it is unlikely that the GMC is undergo- 
ing large-scale, coherent rotation. The position-velocity 
plot for the molecular cloud in Figure 10 (red points) 
has a zero slope, as does NGC 2264, indicating that it 
has no significant linear gradient. It is difficult to say 
whether or not the present velocity field originated dur- 



ing the cloud's formation or during some later stage in its 
evolution. If the GMC originally had a more organized 
velocity field and did not inherent its present field dur- 
ing formation , perhaps a series of interactions with ex- 
ternal forces or internal events — such as turbulence and 
star formation activity — have washed out any signature 
of a systematic velocity gradient which was previously 
present in the molecular cloud. 

The Rosette 

The separation between the gradient directions in the 
Rosette molecular cloud and the surrounding atomic gas 
is approximately 130°. Figure 11 shows that in the veloc- 
ities of the molecular cloud tend to increase from roughly 
7-14 km s~^ in the North-East to as high as ~ 17 
km s~^ in the South and South- West. The gradient in 
the H I is directed perpendicular to the Galactic plane. 
We note that, of the clouds in our sample, the H I maps 
near the Rosette probably suffer the greatest degree of 
line-of-sight blending. The cloud is very close to the 
Galactic plane, and at the distance of the Rosette, 1600 
pc, we have the lowest spatial resolution (^ 17 pc) in the 
atomic gas. Therefore the surface density and velocity 
maps displayed in Figures 5 and 11 almost certainly fail 
to capture many of the local, small-scale variations in the 
structure of the atomic gas. 

We measure a gradient in the molecular cloud of 0.09 
km s~^ pc~^, consistent with the 0.08 km s^^ pc~^ mea- 
sured by Williams et al. (1995). If the gradient is due 
to rotation, the specific angular momentum of the cloud 
is approximately 26 pc km s~^. Based on the estimate 
of the minimum value of the accumulation radius, the 
specific angular momentum in the surrounding ISM is 
nearly 3 times larger. 

Figure 5 shows the outline of the molecular cloud over- 
laid on a surface density map of the atomic gas, which is 
integrated over the range v = 4 27 km s^^. A "shell" 
of H I containing two high-density peaks is associated 
with the molecular cloud. The brightest peak of the shell 
sits at the southern edge of the molecular cloud near a 
latitude of 6 = —3°, and another peak in the surface 
density occurs near b = —1.5°. The southern portion 
of the H I shell appears to mimic the ring-like structure 
in the molecular cloud sitting just to the north of it. 
When the location of the structures in both the molecu- 
lar and atomic gas are compared with the first-moment 
map in Figiue 11, we sec that the southern portion of 
the H I shell is moving near the same velocity, ^ 13 - 
14 km s^^, as the south-east segment of the ring in the 
GMC. The western half of the ring in the GMC is moving 
at higher velocities. Also note how the H I peak locate d 
near b = —1.5° is moving slightly faster 14 km s""'^) 
than the molecular gas in that region 13 km s"'^). All 
of this suggests that both the H I shell and the ring-like 
structure in the GMC are expanding. Kuchar & Bania 
(1993) demonstrated that the H II region NGC 2244 (the 
Rosette Nebula) could have given rise to the expansion 
in the atomic gas. And Williams et al. (1995) found 
that certain properties of clumps within the GMC vary 
with distance from the nebula, centered at / = 206.25°, 
b= -2.11 (Celnik 1983). 

Taken all together, the evidence leads us to suggest 
that the gradient in the Rosette molecular cloud is not 
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caused by rotation, hut by the high-luminosity H II re- 
gion, NGC 2244. High-energy winds of stars in NGC 
2244 may have excavated a hole in the GMC, causing 
H I to be swept up into a high-density ridge by the ex- 
panding molecular gas (e.g., Kuchar et al. 1991, Kuchar 
& Bania 1993). In this picture, because the nebula has 
had less of an impact on distant regions of the cloud, 
these distant regions are moving at lower velocities com- 
pared to molecular gas near the H II region. 

5. IMPLICATIONS FOR GMC FORMATION 

Our key findings are that the angular momentum in the 
GMCs is less than that in the surrounding atomic gas, 
and the velocity gradient position angles in the molecular 
and atomic gas arc widely divergent — with Orion A being 
the one exception. This leads us to suggest that rotation 
may not be the best explanation of the velocity fields 
observed in the GMCs. 

Traditionally, at least three possible solutions have 
been invoked to resolve the angular momentum problem: 

1. One or more of the assumptions in the theory are in- 
appropriate. For instance, if the average surface density 
of the precursor atomic gas is underestimated, this will 
lead to an overestimate in the accumulation radius and, 
consequently, an overestimate of the angular momentum 
initially imparted to a GMC. In our analysis, we used 10 
Mq pc~^ for the mean surface density of the precursor 
gas, twice the average value in the Solar vicinity. Assum- 
ing that the initial gas had an even higher density, say 20 
Mq pc~^, this would not change the main result, namely 
that jHi is consistently greater than j'gmc- This is be- 
cause the accumulation radius depends weakly on the 

initial surface density of the gas {Ra oc Sjjj ). On the 
other hand, if collapsing molecular clouds do not gather 
material far from the Galactic plane, using larger values 
of Ra might be appropriate. Yet this would only exac- 
erbate the discrepancy between the predicted and obser 
ved angular momenta. It might also be argued that using 
the effective radius of a GMC leads to underestimates in 
jco if the GMC is filamentary. In the previous section 
we took this into consideration with Orion A and re- 
calculated jco assuming a cylindrical morphology. This 
raised jcOj which, nevertheless, remained less than jhi- 
Moreover, varying the size of the GMC or the region from 
which it gathers material does not solve the problem of 
the gradient directions in the molecular and atomic gas 
being unaligned. 

2. An alternative explanation of the angular momen- 
tum problem is that though the assumptions regarding 
the initial conditions may be valid, there may be some 
kind of external braking force which rapidly reduces the 
angular momentum of a GMC during its initial condensa- 
tion. Magnetic braking is often evoked as an angular mo- 
mentum shedding mechanism (e.g., Mouschovias 1977; 
Fleck & Clark 1981; Mestel & Paris 1984; Rosolowsky et 
al. 2003). Zeeman splitting of the OH 18-cm line and the 
21-cm line of neutral hydrogen have been used to mea- 
sure the magnetic field strengths of GMCs. Since GMCs 
are magnetized and MHD effects are expected to play a 
significant role in their evolution, magnetic braking -in 
which magnetic field lines anchoring a GMC to the am- 
bient ISM provide the tension necessary to slow down 
rotation — is a possible solution to the angular momen- 
tum problem. Heiles & Troland (2005) measured the 



mean magnetic field strength in the cold neutral medium 
of the Milky Way to be Bo = 6 ± 1.8 ^G. The braking 
time is set by the time it takes for Alfven waves to travel 
across a region of gas having a moment of inertia com- 
parable to the GMC (Mestel & Paris 1984). In order for 
magnetic braking to be efficient at slowing down cloud 
rotation, it would have to occur on timescales no greater 
than the timescale for cloud collapse. For gas having an 
initial density of 1 cm"^ (=0.025 Mq pc~^), the Alfven 
speed is Bo(47rp)~^/^ w 9.5 km s~^. In a region having 
an accumulation radius of 70 pc, this corresponds to a 
braking time of roughly 7.4 Myr. By comparison, the 
timescale for self-gravitational cloud formation in a re- 
gion having the same density, assiuning that this process 
occurs on a timescale close to the dynamical free-fall time 
(Mouschovias & Paleologou 1979; Mestel & Paris 1984; 
Elmegreen 2007), is {Stt/AGpY/'^ « 44 Myr. 

However, the effectiveness of magnetic braking largely 
depends on how the mass of the forming cloud compares 
to its magnetic critical mass. Mestel & Paris (1984) argue 
that braking will efficiently slow down rotation only if the 
mass, M , of the forming cloud is much less than its mag- 
netic critical mass, Mc- Using Crutchcr's (1999) mag- 
netic field strength measurements of molecular clouds, 
McKee & Ostriker (2007) infer that GMCs have approx- 
imately M > 2Mc, that is, they are magnetically su- 
percritical. Elmegreen (2007) argues that the ISM out 
of which GMCs form has a magnetic field that is "near- 
critical," i.e., M ^ Mc- Furthermore, magnetic braking 
slows rotation most efficiently when the cloud's angular 
momentum vector is perpendicular to the magnetic field 
(Mouschovias & Paleologou 1979). Due to the complex- 
ity of making precise measurements of the magnetic field 
in the ISM, it is not established whether or not this is 
the case. Clearly, more observatio ns are needed in order 
to determine the importance of magnetic fields in GMC 
evolution. 

3. Finally, perhaps GMCs are not rotating after all, 

and the "problem" is unfounded. As some authors have 
pointed out, the interpretation that velocity gradients 
indicate rotation is not unique. Expansion and shear, 
for instance, also produce velocity gradients. The nu- 
merical simulations of Hennebelle et al. (2008) suggest 
that the converging flows of atomic gas could produce 
GMCs. If shearing occurs at the interface of colliding 
flows, might this resiilt in an excess of shear, that is, 
values of O that are higher than the shear arising from 
the local Galactic rotation curve? Assuming a flat ro- 
tation curve, the local value of in the solar neighbor- 
hood is 0.025 km s~^ pc~^. This is slightly less than, 
but basically comparable to the gradient magnitudes we 
measure in Perseus, Orion A, and MonR2 (Table 3). As 
discussed previously, we arc; underestimating in every 
case, since we have no information regarding the inclina- 
tion of clouds. Also, since bl ending of emission along the 
line-of-sight is most extreme in the cases of NGC 2264 
and the Rosette, which are situated close to the Galactic 
plane, it is likely that fifji is underestimated to an even 
greater degree in these cases. Nevertheless, we cannot 
conclude whether or not shear is the cause of the veloc- 
ity gradients based on these observations. It is an issue 
we would like to further investigate with a larger sample 
of clouds. 
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Burhert & Bodcnhcimcr (2000) show that turbulence 
may also cause linear velocity gradients. They demon- 
strated that the gradient magnitude of turbulent cores 
scales with size as f2 oc Figure 14 shows ft 

as a function of R for both the molecular clouds and 
H I. Neither the GMCs nor the H I appear to follow the 
O oc R^'^-^ relationship, though we do not have enough 
data points to make a definite conslusion. 

Studying the properties of GMC kinematics in other 
galaxies will also help to develop our picture of GMC 
formation. In their analysis of 45 GMCs in M33, 
Rosolowsky et al. (2003) did the first systematic, ex- 
tragalactic study of GMC angular momentum proper- 
ties. They showed that simple GMC formation theories 
consistently overestimate the magnitude of the observed 
angular momentum. They measured the velocity gradi- 
ents from high-resolution ^^C0( J = — >■ 1) data, finding 
that the gradients of M33 clouds are similar in magnitude 
to Galactic clouds. They then tested several formation 
models by calculating the accumulation radii of the cat- 
aloged clouds which are predicted by the various models, 
including the Toomre and Parker instabilities. On av- 
erage, the theories, which do not include the effects of 
magnetic fields on rotation, overpredict observed magni- 
tudes of velocity gradients by more than a factor of 5. 
And 40% of the GMCs in M33 are counterrotating with 
respe ct to the sense of galactic rotation. 

We extend this study by performing the analysis estab- 
lished here on a much larger sample of GMCs in M33 cat- 
aloged by Rosolowsky et al. (2003). Since M33 has a rel- 
atively low inclination of ~ 51° (Corbelli & Salucci 2000), 
such a study has the advantage of bypassing the prob- 
lem of source confusion along lines of sight that arises 
when doing comparable surveys of molecular clouds in 
the Milky Way. 

6. SUMMARY AND CONCLUSIONS 

We have presented a detailed comparison between the 
kinematics of five Galactic molecular clouds and the 
atomic gas that surrounds them. We developed a method 

for selecting regions of H I that are associated with the 
GMCs and found that each GMC was nearby high- 
density H I which peaked above the mean Galactic value. 
First-moment maps were created using the ^"^CO emis- 
sion of the GMCs and the H I 21-cm emission, and then 
a plane was fit to each map of velocity centroids. We de- 
termined the magnitudes and directions of the velocity 
gradients from the coefficients to the fits. From these ob- 
servations and measurements, we arrived at the following 
conclusions: 

1. Orion A, Perseus, and the Rosette each have a sig- 
nificant linear velocity gradient across the face of the 

molecular clouds, while NGC 2264 and MonR2 have com- 
plex, non-linear velocity fields. The Perseus molecular 
cloud has the strongest linear gradient in the sample, 
with a magnitude of 0.23 km s^^ pc"-"^. 

2. The atomic gas associated with Orion A, MonR2 
and the Rosette has significant linear velocity gradients, 
regardless of whether the molecular cloud has one. The 
HI gradients range from 0.019 to 0.067 km s"'^ pc^^, 
or 0.76 to 2.7 times the shear in the solar vicinity as 
measured by the Galactic rotation curve. 

3. If the gradients in the molecular and atomic gas 
were due to rotation, the specific angular momentum in 



the GMCs is less than that predicted by the formation 
scenario in which a GMC preserves angular momentum 
while undergoing top-down collapse from the surround- 
ing ISM by a factor of 1 to 6. The discrepancy can be 
narrowed if different assumptions are made regarding the 
initial density and geometry of the gas. But the direc- 
tion of the trend that the observed angular momentum 
is less than the predicted — remains the same. 

4. We observe large differences between the velocity 
gradient directions in the GMCs and the atomic gas, with 
Orion A being the one exception. At more than 130°, 
two of the most extreme angle separations can be seen 
in the MonR2 and the Rosette systems (Figm-e 10 and 
11). Fiirthermore, the gradient directions in neither the 
molecular nor the atomic gas are in alignment with the 
overall direction of Galactic rotation. If the gradients 
were due to rotation, this indicates that some GMCs are 
counterrotating with respect to the Galaxy. 

5. That the velocity gradient position angles in the 
atomic and molecular gas are divergent indicates that 
the GMCs in our sample probably did not inherit their 
present velocity fields from the atomic gas from which 
they formed. Finally, in at least two cases, NGC 2264 
and the Rosette, a good explanation of the morphol- 
ogy and kinematics observed in the gas is that they are 
caused by stellar winds from O stars in the H II regions 
located in the GMCs, not rotation. 
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TABLE 1 
GMC Properties 





Coordinates 


Distance 


Mi3 


MqmC 


Projected Area 




v{H I) 


References 


Cloud 


^0 


bo 


[pc] 


[10* Mq] 


[104 Mq] 


[pc2] 


[km s^^] 


[km s^^] 


Distance ^^CO 


Perseus 


159 


-20 


320 


1.4 


1.2 


825 


7.5 


-8,14 


1 6 


Orion A 


210 


-19 


414 


2.4 


6.9 


1510 


9.6 


-7, 22 


2 7 


NGC 2264... 


202 


1 


800 


2.5 


2.2 


674 


6.6 


-7, 15 


3 8 


MonR2 


213 


-13 


830 


3.0 


9.0 


2940 


10.8 


-5,23 


4 9 


Rosette 


207 


-2 


1600 


4.0 


7.7 


2200 


13.1 


4,27 


5 10 



TABLE 1 

Properties are determined using the methods outlined in §3 of the text. Masses calculated using -'-■^CO observations are 

LISTED under M13. THE GMC MASSES CALCULATED BY PREVIOUS AUTHORS ARE LISTED UNDER MqmC- 

References.-(I) de Zeeuw et al. 1999; (2) Genzel et al. 1981; (3) Sagar & .losiii ; (4) Racine 1968: (5) Blitz & Stark 1986; 
(6) Sargent 1979; (7) Menten et al. 2007; (8) Blitz 1978; (9) Maddalena et al. 1986; (10) Williams et al. 1995. 



TABLE 2 

Accumulation Radii 



Cloud 




Ra 


Perseus 


16.2 


42.4 


Orion A 


21.9 


76.0 


NGC 2264... 


14.6 


51.9 


MonR2 


30.6 


83.6 


Rosette 


26.5 


78.8 



TABLE 2 

The effective GMC radii, -RqmCi and the accumulation radii, R^, are listed in units of pc. Lower limits of _Ra are 
estimated using Equation 12, assuming that the GMCs initially formed from H ihaving an average surface density of 10 

Mq PC~2. 



TABLE 3 
GMC Properties: Dynamics 

Cloud f^GMC S^Hi 6'gmc ''ni icivic Jhi 

[0.01 km s~^ pc~^] [0.01 km s~^ pc~^] [deg] [deg] [pc km s'^] [pc km s'^] 

Perseus 23.1 6.70 -20.2 99.4 24.2 60.3 

Orion A 22.0 3.71 -111 -102 42.2 [84.1] 107 

NGC 2264... 4.60 1.87 47.8 -173 3.92 25.2 

MonR2 6.68 3.48 -72.4 -11.2 25.5 120 

Rosette 9.18 2.36 -132 1.12 25.8 73.4 



TABLE 3 

The magnitude of the velocity gradient, H, and the direction of rotation measured in degrees East of North are 
determined using the methods outlined in the text. The specific angular momentum, j, IS calculated for each GMC and 

H IREGION assuming that the GRADIENTS ARE DUE TO ROTATION. LISTED IN BRACKETS IS j FOR OrION A, ASSUMING THAT THE GMC 
has THE MORPHOLOGY OF A CYLINDER WITH A ROTATION AXIS PERPENDICULAR TO ITS LONG AXIS. TYPICAL ERRORS ARE 

SCI = 0.001 - 0.005, 56»(CO) = 0.5 - 2°, 5e»(HI) = 5-10°, and Sj =0.5-2 PC km S"!. 
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Fig. 1. — Orion A. The bottom figure plots the average spectra of the '^^CO emission (sohd line) in the GMC and of the H I emission 
(dashed line) in the region of atomic gas. The top figure shows a surface density map of H I overlaid with an outline of the molecular 
cloud. The range of H I surface density, in units of Mq pc~-^, is marked in the top left-hand corner, and the contour spacing is 1 Mq pc~^. 
Atomic gas within the dashed line is used for the subsequent analysis described in the text. 
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Fig. 2. — Perseus. Same as Figure 1. 
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NGC 2264. Same as Figure 1. 
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Fig. 4. — MonR2. Same as Figure 1. 
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Fig. 5.— 



The Rosette. Same as Figure 1. 
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Fig. 6. — Orion A. The top figure shows the intensity-weighted first moment maps of H I (grayscale) overlaid with ^^CO (color). Both 
maps have contour spacings of 0.5 km s~^, and the velocity ranges of the maps are indicated by the color bars. The first moment map 
to the right shows the H I without the ^^CO overlaid. The bottom figure plots the central velocity versus perpendicular offset from the 
rotation axis for pixels in the above map located within the dashed line. The ^^CO data (in red) are binned every 3 pc and the H I data, 
every 6 pc. The error bars indicate the dispersion within the bins. 
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Fig. 7. — Perseus. Same as Figure 6, except the H I data in the position-velocity plot are binned every 4 pc. The non-linearity of the 
H I position-velocity plot indicates that there is not a linear velocity gradient over the entire region within the dashed circle. In the following 
figure, we see that by changing the reference position, there is a significant linear gradient in the atomic gas. 
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Fig. 8. — Perseus. Same as Figure 6, accept the H Idata in the position-velocity plot are binned every 4 pc. The surface density map 
in Figure 2 shows that a high-density H Ipeak is located to the West of the molecular cloud. By changing the reference position to 
Iq, bo = 157.5, —18, near the center of the peak, we find that there is a significant linear gradient across the field centered on the peak. 
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Fig. 9. — NGC 2264. Same as Figure 6, except the H I data in the position-velocity plot are binned every 9 pc. Although there is no 
significant linear gradient across the face of the GMC, we plot the "rotation axis" (red line) measured using Equation 10 as if there was. 
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Fig. 10. — MonR2. Same as Figure 6, except the H I data in the position-velocity plot are binned every 9 pc. Although there is no 
significant linear gradient across the face of the GMC, we plot the "rotation axis" (red line) measured using Equation 10 as if there was. 
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Fig. 11. — Rosette. Same as Figure 6, except the H I data in the position- velocity plot are binned every 17 pc. 
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Fig. 12. — Specific angular momentum of the atomic hydrogen, jhi, versus GMC specific angular momentum, j'gmc- Overplotted is a 
least-squares fit to the data: ;?hi oc JG^fc^"^"' 
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Fig. 13. — The gradient position angles of the atomic hydrogen, i/ihi, versus position angles of the GMCs, Vgmc appear to be uncorrelated. 
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Fig. 14. — The gradient magnitude versus size for H I and GMCs (red points). Ovcrplottcd is the proportionality f! oc i? "'^ that Burkcrt 
k Bodenheimer (2000) found for turbulent cores. 
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